Double Neutron Star Binaries: Implications for LIGO 
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Double neutrons are especially important because they give most accurate informations on the 
masses of neutron stars. Observations on double neutron stars show that all masses of the neutron 
stars are below l.5M@. Furthermore, two neutron stars in a given double pulsar are nearly equal 
in mass. With hypercritical accretion, we found that the probability of having companion mass 
> 1.5M© is larger than 90%, while there is no observations on such systems. We believe that those 
companions with masses higher than 1.5-Mq went into black holes, which is consistent with our 
preferred maximum neutron star mass M^g x sa 1.5Mq due to the kaon condensation. In this work, 
we point out that the black-hole neutron star binaries are 10 times more dominant than double 
neutron star binaries. As a result, black-hole, neutron star binaries can increase the LIGO detection 
rate by a factor 20. 

PACS numbers: 97.60.Bw, 97.60.Lf 



I. INTRODUCTION 



The observational evidence as well as the calculational 
basis, is strong that neutron star (NS) binaries evolve 
from double helium stars, avoiding the common envelope 
evolution of the standard scenario [1|. (A helium star 
results in a giant when the hydrogen envelope is lifted 
off - in binary evolution by being transferred to the less 
massive giant in the binary.) In the standard scenario of 
binary NS formation after the more massive giant trans- 
fers its hydrogen envelope to the companion giant, the 
remaining helium star burns and then explodes into a 
NS. In about half the cases the binary is not disrupted 
in the explosion. The NS waits until the remaining gi- 
ant evolves (and expands) in red giant following its main 
sequence hydrogen burning. Once the envelope is close 
enough to the NS, the latter couples to it hydrodynami- 
cally through gravity. Some of the material in the enve- 
lope is accreted onto the NS, although most flies by in the 
wake, being heated in the process, and is lost into space. 
The energy to expel the matter comes from the drop in 
potential energy as the orbit of the NS tightens. Formu- 
las for the tightening and the amount of mass accreted 
by the NS were given by Bethe & Brown 0. 

Chevalier first estimated that in common envelope evo- 
lution the NS would accrete sufficient matter to evolve 
into a black hole (BH) Brown suggested the double 
helium star scenario, in order to save the first born NS 
P|. In this scenario, mass exchange of the hydrogen en- 
velope takes place while both stars burn helium. There 
is not sufficient time for this mass to be accepted by ei- 
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ther star p| and it is lost into space, leaving a binary 
of helium stars. The two giant progenitors must have 
main sequence masses within ~ 4% of each other, in or- 
der to burn helium at the same time, a highly restrictive 
requirement. However, the observations of nearly equal 
masses of the two NS's within the binaries give support 
to this scenario, as we show. 



The above scenario for binary evolution was made 
quantitative by Bethe & Brown [2( who calculated that 
in a typical case, the NS would accrete ~ 1M , taking it 
into an ~ 2.4M Q low mass black hole (LMBH), similar 
to the BH we believe resulted from SN1987A, although 
somewhat more massive than the latter. Thus, if the 
NS had to go through common envelope evolution in a 
hydrogen envelope of ^> 1OM from the giant, it would 
accrete sufficient matter to go into a LMBH. Therefore, 
when the giant evolved into a helium star, which later 
exploded into a NS, a LMBH-NS binary would result 
provided the system was not broken up in the explosion. 
(About 50% of the time the system survives the explo- 
sion.) In this work, we extend Bethe & Brown work to 
include the hypercritical accretion during both red-giant 
and super-giant stages of the second star which evolves 
later. 



The above scenario was estimated to take place 10 
times more frequently than binary NS formation which 
required the two stars to burn helium at the same time, 
and, because of the greater mass of the BH, the mergings 
of binaries with LMBH to be twice as likely to be seen 
as those with only NS's. This is the origin of the factor 
20 enhancement of gravitational mergers to be observed 
at LIGO, over the number from binary NS's alone. 



II. DOUBLE NEUTRON STARS 



Case 1: AT > 10%, P ~ 90%, AM ~ 0.9M P 



We list in Table [I] the 5 observed NS binaries with 
measured masses. Note that they are consistent with our 
preferred maximum mass of neutron star M™g x = 1.5M Q 
due to the kaon condensation |l2j| . In addition, the two 
NS's in a given binary have very nearly the same mass, 
as would follow from the double helium star scenario, as 
we discuss in next section. 
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TABLE I: 5 observed NS binaries with measured masses. 



Object 



Mass (M ) Companion Mass (M ) Refs. 
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1 907+0. 005 
l.OO/ _ Q 005 


tl orn+OOO 5 








f J0737 
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The very nearly equal masses of pulsar and companion 
in B1534+12 and B2127+11C is remarkable. We show 
below that B1913+16 comes from a region of giant pro- 
genitors in which the masses could easily be as differ- 
ent as they are. The uncertainties in J1518+4904 are 
great enough that the masses could well be equal. The 
the double pulsar J0737-3039A and J0737-3039B were 
probably very nearly the same before a common envelope 
evolution in which the first formed NS J0737— 3039A ac- 
creted matter from the evolving (expanding) helium star 
progenitor of J0737— 3039B (Case 2 in next section) in 
the scenario of Dewi & van den Heuvel [13 . 



III. FATE OF COMMON ENVELOPE 
EVOLUTION 

In the standard scenario the first-born NS would go 
through common envelope evolution with a giant which 
must have ZAMS mass at least ~ IOMq if it is to later 
end up as a NS .1]. In Fig. three typical cases of bi- 
nary NS evolution are summarized. In these estimates, 
we assumed that both the birth rate and the life time are 
proportional to M^^s- With these assumption, the 4% 
difference in the ZAMS mass corresponds to 10% differ- 
ence in the life time. Hence the population probability 
for the ZAMS masses to be within 4% difference in mass 
is about 10%. 

• Case 1: [90% probability] This corresponds to the 
case with initial mass difference AMzams > 4%, 
with life time difference AT > 10%. The first born 
NS can accrete both in red giant and in super giant 
stage of the second star which evolves later. Due to 
the hypercritical accretion, AM = 0.9M Q (0.7M Q 
in giant stage, 0.2M Q in supergiant stage) can be 
accreted to the first-born NS. The accreted mass 



Case 2: AT < 10%, P ~ 10%, AM ~ O.2M 
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Case 3: AT < 1%, P < 1%, AM ~ 0.0Af G 
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FIG. 1: Three typical cases of binary NS evolution. AT is 
the difference in the life times of binary progenitors, P is 
the formation probability of binary progenitors, AM is the 
accreted mass by the first-born NS. 



AM was estimated using the formula given in the 
Appendix of Belczynski et al. [l4| . 

• Case 2: [10% probability] This corresponds to the 
case with initial mass difference AMzams < 4%, 
but not close enough to burn the helium at the same 
time. AM = 0.2M© can be accreted to the first- 
born NS during the supergiant stage of the second 
star. 

• Case 3: [< 1% probability] This corresponds to 
the case in which the initial mass differences are so 
close to burn the helium at the same time. Nothing 
can be accreted because two NS's are formed almost 
at the same time. 

As in Fig. |3 if the NS's are to be born with initial 
masses between 1.2M and 1.5M Q as seen in double 
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delay. Although the observational indication of high NS 
mass is strong, this mass would be brought down with 
the 4/3 power of white dwarf mass if the latter were in- 
creased, and still fit the same period change. Thus, we 
believe the case for such a massive NS to only be settled 
with a sufficiently accurate measurement of the Shapiro 
decay which pins down the white dwarf mass. It should 
be noted that just in the evolution of NS, white-dwarf 
binaries there is ample possibility for substantial accre- 
tion from the evolving progenitor of the white-dwarf, so 
these binaries are the place to look if one wants to find a 
high-mass NS. 



IV. OBERVABILITY PREMIUM FOR LMBH-NS 
BINARIES 



FIG. 2: Companion (first-born NS) mass in double NS's and 
their formation probability. The area within the square indi- 
cates the probability. Filled circles are observed companion 
masses in double NS's. Here, NS's are assumed to be born 
with initial masses between 1.2Mq and 1.5Mq. Note that 
the probability of having higher mass (> 1.5Mq) companion 
is more than 90%. We believe that they went into black holes 
and the maximum mass of NS M^g x — 1.5Mq is consistent 
with observations. 



NS's, the first-born NS's in Case 1 result in 2.1 - 2.4M Q 
due to the hypercritical accretion. Similarly, those in 
Case 2 result in 1.4 — 1.7M Q . Furthermore, one shouldn't 
forget that these massive pulsars in Case 1 would be co- 
pious, because there is no special condition of the giant 
progenitors of the binary NS having to be within 4% of 
each other in ZAMS mass. In fact, the binaries with 
pulsar mass 2.1 — 2.4Af would be a factor of 10 more 
frequent than those with companion masses 1.2 — 1.7Mq. 
Since we don't see pulsars with such high masses, we be- 
lieve that they must have gone into LMBH's. 

Our argument in this note alone does not exclude NS 
in the mass range up to 2.1M Q , but as mentioned ear- 
lier the Bethe & Brown |l6j | argument that the maximum 
mass of the NS in 1987A, which we believe went into a 
LMBH, of 1.57M Q further constrain the maximum NS 
mass, if our belief is correct. In Fig. [21 we draw the ex- 
pected maximum mass of NS M™c x = 1.5M Q which was 
estimated from kaon condensation [l^. Note that the 
probability of having higher mass (> 1.5Mq) compan- 
ion is more than 90% while there is no observations on 
such systems. We believe that this indicates that those 
NS's with masses > 1.5M Q went into black holes after 
common envelope evolution. 

The NS mass in the helium white-dwarf, NS binary 
J0751-1807is quoted as 2.1±O.2OM [l|. The NS mass 
in J075 1 — 1807 is measured from the period change due to 
gravitational wave emission. The companion white dwarf 
mass is constrained by a marginal detection of Shapiro 



Why haven't we seen any LMBH-NS binaries? Van 
den Heuvel 0,0] has pointed out that NS's form with 
strong magnetic fields 10 12 to 5 x 10 12 gauss, and spin 
down in a time 

r sd ~ 5 x 10 6 years (1) 

and then disappear into the graveyard of NS's. The 
relativistic binary B1913+16 has a weaker field B ~ 
2.5 x 10 10 gauss and therefore emits less energy in mag- 
netic dipole radiation. Van den Heuvel estimates its spin- 
down time as 10 8 yrs. There is thus a premium in obser- 
vational time for lower magnetic fields, in that the pulsars 
can be seen for longer times. Taam & van den Heuvel 
|19| found empirically that the magnetic field of a pulsar 
dropped roughly linearly with accreted mass. This ac- 
cretion can take place from the companion in any stage 
of the evolution. A pulsar that has undergone accretion 
is said to have been "recycled". In B1913+16 the pulsar 
magnetic field is ~ 2.5 x 10 10 gauss and in B1534+12 
it is - 10 10 gauss. In J0737-3039A it is only 6.3 x 10 9 
gauss. These "recycled pulsars" will be observable for 
> 100 times longer than a "fresh" (unrecycled) pulsar. 

The same holds for LMBH-NS binaries. The NS is 
certainly not recycled, so there is an about 1% chance 
of seeing one as the recycled pulsar in a binary NS. But 
we propose 10 times more of the LMBH's than binary 
NS's, of which we observe 5. Thus the total probability 
of seeing the LMBH binary should be about 50%. How- 
ever, there may be additional reasons that the LMBH-NS 
binary is not observed. 

In Table |Hj the estimated LIGO detection rates are 
summarized. For completeness we added the contribu- 
tion from BH-BH mergers obtained by Portegies Zwart 
& McMillan |20j . They have suggested a large number 
of gravitational mergings of high-mass BH's ejected from 
globular clusters. Their predictions should be tested rel- 
atively early in the LIGO development. 
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TABLE II: Predicted LIGO Detection Rates (yr~ 



Binary Type LIGO I LIGO II Chirp Masses (M ) 



NS-NS 1 " 


0.0348 


187 


1.0 - 


1.3 


BH-NS n 


0.696 


3740 


1.3 - 


2.7 


BH-BH** 


0.58 


2450 


rv 


6 


Total 


1.31 


6377 







t NS-NS detection rates are from Kalogera et al. 21]. ^ 
BH-NS detection rates are obtained by multiplying factor 20 
to NS-NS detection rates from Bethe & Brown [2j. ** 

BH-BH detection rates are from Portegies Zwart & 
McMillan [^J with the modification of BH mass 7M Q . 



CONCLUSION 



The Chirp mass, which will be detected with an esti- 
mated accuracy of ~ 0.002M© 



M, 



chirp 



M 3/5 M 2/5 



(MiM 2 ) 3 / 5 

(Mi + M 2 y/5 



(2) 



would be 1.22M for the merging of two IAM & NS's, 
for the merging of the least massive 2.1M Q BH, result- 
ing from the common envelope evolution of a 1.4M NS 
in a ZAMS 1OM progenitor companion the chirp mass 
would be 1.49Mq. Estimated Chirp masses are summa- 
rized in Table ILTI We may have to wait for LIGO which 
will be able to measure "chirp" masses quite accurately. 
The chirp mass of a NS binary should concentrate near 
1.2Mq, whereas the LMBH-NS systems should have a 
chirp mass of > 1.4M©, and there should be ~ 20 times 
more of the latter. 



The discovery of the double pulsar increased estimated 
rate for gravitational merging by a factor of 6—7 over that 
of Kim et al. [22j. In addition to these effect, we find 
10 times more LMBH-NS binary mergings, with larger 
chirp mass than NS-NS binaries because of the accretion 
in forming the BH so that these mergings multiply the 
binary NS ones by a factor of 20. According to our esti- 
mates, LIGO I would be able to detect one mergings per 
year. 
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